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Abstract. We describe updated calculations of QQ production in pp and 1C p inter- 
\ actions. We compare these results to total cross section data and discuss how the 

C*~) ■ baseline cross sections extrapolate to heavy ion collider energies. We touch upon the 

C*~) | differences between leading and next-to-leading order heavy quark production. Fi- 

nally, we discuss the implications of our calculations for quarkonium production. Our 
discussion here focuses on bottom quarks. 



in 
o 

(N 
O 

a: 



- 1—1 

X 



Keywords: heavy quarks, perturbative QCD 
PACS: 25.75.-q 



a: 

^ ■ 1. Introduction 



Heavy quark production is of great relevance for high energy nuclear collisions. The large 
heavy quark masses means that their production can be calculated in perturbative QCD. 
They are produced in the initial nucleon-nucleon collisions and, as such, provide valuable 
information about the early state of the system. It is therefore important to make systematic 
studies to obtain the fullest possible information. The baseline rates in pp interactions are 
essential to determine the expected total cross sections and the unmodified quark distribu- 
tions. 

There are many effects which can be studied through systematic heavy quark measure- 
ments. Heavy quark decays are expected to dominate the lepton pair continuum from the 
7/\|/(cc) and J(bb) up to the mass of the Z° [ U §, [|. Thus the Drell-Yan yield and any 
thermal dilepton production will essentially be hidden by the heavy quark decay contribu- 
tions [ [j]]. The shape of the charm and bottom contributions to this continuum could be 
significantly altered by heavy quark energy loss Since B mesons decay to the 

secondary J/\\t spectrum could be modified by energy loss effects on the primary B [ 2[. If 
the heavy quark energy loss is large, it may be possible to extract a thermal dilepton yield 
if it cannot be determined by other means [ 6fl. Heavy quark production in a quark-gluon 
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plasma has also been predicted [ pj. This additional yield can only be determined if the 
relative AA/pp rate can be accurately measured. Finally, the cc yield is a useful reference 
for J/\\i production since a (//ty) / (cc) enhancement has been predicted due to cc recom- 
bination [f| | g, This recombination requires more than one cc pair produced in an 
event, easily reached in heavy ion collisions at RHIC and LHC. 

An update to the pp baseline is relevant now because the gluon distributions have 
changed significantly since the calculations of Ref. [ [Tl| ], The rise at low momentum frac- 
tions, x, is not as large as previously expected, significantly lowering the cc cross sections 
in particular. Less attention has been paid to the bb cross sections, perhaps due to a lack of 
pp data. However, since a new measurement has recently become available [ p2[ ], we will 
try to make up for this shortcoming here by concentrating on bottom production. 



2. QQ Total Cross Sections 

At leading order (LO) heavy quarks are produced by gg fusion and qq annihilation while 
at next-to-leading order (NLO) qg and qg scattering is also included. To any order, the 
partonic cross section may be expressed in terms of dimensionless scaling functions ft- ' ' 
that depend only on the variable r| [ pj| , 

m k=0 1=0 \ m Q ) 

where s is the partonic center of mass energy squared, mg is the heavy quark mass, p is the 
scale and r| = s/Am 2 Q — 1. The cross section is calculated as an expansion in powers of (X s 
with k = corresponding to the Born cross section at order o(a^). The first correction, 
k=l, corresponds to the NLO cross section at o(oc,). It is only at this order and above 
that the dependence on renormalization scale, pr, enters the calculation since when k = 1 
and I — 1, the logarithm ln(// 2 /mg) appears. The dependence on the factorization scale, 
fjf, the argument of tt s , appears already at LO. We assume that /jr = fip = /j. The next-to- 
next-to-leading order (NNLO) corrections to next-to-next-to-leading logarithm have been 
calculated near threshold [ |l3|| but the complete calculation only exists to NLO. 

The total hadronic cross section is obtained by convoluting the total partonic cross 
section with the parton distribution functions (PDFs) of the initial hadrons, 




where the sums i and j are over all massless partons and x\ and are fractional momenta. 
The PDFs, denoted by FP, are evaluated at scale fj. All our calculations are fully NLO, 
applying NLO parton distribution functions and the two-loop evaluation of a s to both the 
(a 2 ) and O (ocf ) contributions, as is typically done n4Q. 

To obtain the pp cross sections at the RHIC and LHC proton and ion energies, we first 
compare the NLO cross sections to the available cc and bb production data by varying the 
mass, niQ, and scale, p, to obtain the 'best' agreement with the data for several combinations 
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of mg, /j, and PDF. These 'best' fits aim to describe the available data without a scaling 
factor, the experimental 'K' factor. We use the recent MRST HO central gluon [ pa], 
CTEQ 5M [ and GRV 98 HO [ fj] distributions. The results for the bb cross section 
in %~p interactions, where there are the most data, are shown in Fig. |l]for the MRST HO 
densities with the SMRS pion densities [ n8[j. We calculate the NLO cross sections for 
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Fig. 1. Total bb cross sections in K~p interactions at fixed-target energies as a function of 
the bottom quark mass. All calculations are fully NLO using the MRST HO (central gluon) 
parton densities. From the upper left, the plots show results with the renormalization and 
factorization scales equal to mi,/ 2, ntf, and 2mi, respectively. In each plot, from top to 
bottom the curves are nib = 4.25, 4.5, 4.75, and 5 GeV. The bb data can be found in Ref. [ 




4.25 < nib < 5 GeV with scale choices clockwise from the upper left: p = Mb/2, m.b, ar >d 
2nib- The cross sections decrease as p increases because a s (m/,/2) > a 5 (m&) > a s (2m/,) by 
the running of a s . Evolution of the PDFs with p tends to go in the opposite direction e.g. 
the sea quark and gluon distributions rise as x decreases while p increases. At higher scales 
the two effects tend to compensate and reduce the scale dependence but the charm quark 
mass is not large enough for this to occur. The results for bottom quarks are somewhat 
better since the mass is larger. 

We find reasonable agreement with all three PDFs for m/, = p = 4.75 GeV, = /i/2 = 
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4.5 GeV (dashed), and m/, = 2p = 5 GeV, shown in the right hand side of Fig. ^| The 
MRST results cluster together and lie a bit higher than the GRV 98 results for K~p pro- 
duction while the opposite is true for pp production. The three data points from pp — ► bb 
measurements lie somewhat closer together although the E789 point (the square on the left 
hand side of Fig. |2|) lies considerably below the other two points. In fact, it agrees best with 
the NLO calculations using a 'standard' b quark mass of 4.75 GeV. 
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Fig. 2. Total bb cross sections in pp and K~p interactions compared to data. All calcula- 
tions are fully NLO. The curves are: MRST HO (central gluon) with m b = p. = 4.75 GeV 
(solid), irib = p/2 = 4.5 GeV (dashed), and m& = 2p — 5 GeV (dot-dashed); as well as the 
GRV 98 HO distributions with m h = p = 4.75 GeV (dotted), m b = p/2 = 4.5 GeV (dot-dot- 
dot-dashed), and mi, = 2p = 5 GeV (dot-dash-dash-dashed). The pp data are from the E789 
(square) [113], E771 (diamond) [ ^(J and preliminary HERA-B (circle) [ Q collaborations. 

The somewhat larger spread in the n~p calculations may be because the 7t~ PDFs are 
not very well known. The last evaluations, SMRS [ |8[], Owens-7t [ [HJl, and GRV-Tt [ |2| 
were 10-15 years ago and do not reflect any of the latest information on the low x behavior 
of the proton PDFs, e.g. the distributions are all flat as x — > with no low x rise. The GRV 
distributions are based on dynamical parton densities with low initial scale which generate 
larger densities at low x while depleting them at high x, reducing the bb cross sections 
at low y/s relative to the flat distributions of SMRS, see Fig. ||. These pion evaluations 
also depend on the behavior of the proton PDFs used in the original fit, including the 
value of Aqcd- Thus the pion and proton PDFs are generally incompatible. The typical 
x values of bb production are large, 0.14 < x — 2p/^/s < 0.9. At the fixed target energies 
of bb production, qq annihilation dominates while gg fusion is still most important for 
cc production [ \23f\ . The valence-valence u K u p contribution dominates since the valence 
distributions are greatest at large x. 

For charm production, the 'best' agreement with p. = m c is for m c — 1 .4 GeV and m c = 
1 .2 GeV is the best choice for p = 2m c for the MRST HO and CTEQ 5M distributions. The 
best agreement with GRV 98 HO is p — m c = 1.3 GeV while the results with p = 2m c lies 
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below the data for all m c . All five results agree very well with each other for pp^ cc. There 
is also more of a spread in the K~ p — ► cc results. The %~ p — > cc cross sections are a bit 
lower than the data compared to the pp cross sections, suggesting that lighter quark masses 
would tend to be favored. The reason is because the low x rise in the proton PDFs depletes 
the gluon density for x > 0.02 relative to a constant at x — ► for fx = /jq the initial scale 
of the PDF. The K~ p data are also in a relatively large x region, 0.1 < x = 2fx/y/s < 0.3, 
where this difference is important. For the cc results, see Ref. [|24|. 

Before calculating the QQ cross sections at nuclear colliders, some comments need to 
be made about the validity of the procedure. Since the cc calculations can only be made 
to agree with the data when somewhat lower than average quark masses are used and even 
the pp — > bb data suggest ntj, should be smaller, it is reasonable to expect that higher order 
corrections beyond NLO could be large. Indeed, the preliminary HERA-B cross section 
agrees with the NNLO-NNLL cross section in Ref. [ Km, suggesting that the next order 
correction could be nearly a factor of two. Thus the NNLO correction could be nearly as 
large as the NLO cross section. 

Unfortunately, the NNLO-NNLL calculation is not a full result and is valid only near 
threshold. The p~p data at higher energies, while not total cross sections, also show a large 
discrepancy between the perturbative NLO result and the data, nearly a factor of three 
[ ^5]]. This difference could be accounted for using unintegrated parton densities [ 
although these unintegrated distributions vary widely at this point. The problem is then 
how to connect the regimes where near-threshold corrections are applicable and where 
high-energy, small x physics dominates. The problem is increased for charm where, even 
at low energies, we are far away from threshold. 

Our method is perhaps the easiest thing to try-using NLO only and ignoring higher- 
order corrections to fit the data. This is not difficult for cc because the data are extensive. 
However, there is less bb data to go on. The p — ► cc data tends to favor lighter masses. It 
is difficult to say if the same is true for bb. A value of mi, = 4.75 GeV, which underpredicts 
the Tevatron results compared to NLO cross sections [ |25[ ], agrees reasonably with the 
average of the n~ p data. However, for the HERA-B measurement to be compatible with a 
NLO evaluation, the b quark mass would have to be reduced to 4.25 GeV, a value which 
might be more compatible with the Tevatron results. Therefore, if the NNLO cross section 
could be calculated fully, the results would likely be more compatible with a larger quark 
mass. A quantitative statement is not possible at this time. 

If we then assume that the NNLO and higher orders to not have a substantially different 
energy dependence than the LO and NLO results, then we will be in the right ballpark 
at collider energies. If the LO and NLO matrix elements are both evaluated with NLO 
PDFs and the two-loop oc s , the theoretical K factors have a relatively weak s/s dependence, 
< 50% for 20 GeV < < 14TeV. The produced heavy quark distributions might be 
slightly affected since the shapes are somewhat sensitive to the quark mass but charm is far 
enough above threshold at ion colliders for the effect to be small. A somewhat larger effect 
might be expected for bottom. 

We now extrapolate our results to RHIC and LHC energies. The bb cross section is 
shown in Fig. [|. The spread in the bb cross sections is small, ~ 20 — 30% at the ion collider 
energies. The cc cross sections, on the other hand, differ by a factor of two at 5.5 TeV. 
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Defining the theoretical K factor as NLO/LO, where both the numerator and denomi- 
nator are calculated with NLO PDFs and the two-loop a s , is not the only way to proceed. 
It is, however, the most common because, in general, one wants to determine the size of the 
corrections. In principle, it is most correct to use LO PDFs and the one-loop a s with the 
LO matrix elements and NLO PDFs and the two-loop a s only with the (aj) corrections. 
Since the QCD scale, A, is tuned to fit all available data in global analyses of parton distri- 
butions, whether at LO or NLO, the number of loops in the evaluation of oc s is important. 
When a fully LO calculation of the charm cross section is done with a one loop evaluation 
of a s in the global analyses of the parton densities, as with the MRST LO set, the LO cross 
sections increase by ~ 60%. This change is almost solely due to the difference between 
the one and two loop evaluations of oc s . The MRST LO A is 0.204 GeV when rif = 3, 
leading to a] loop = 0.287 and a? loops = 0.220 for /j = 2m c . The A associated with the 
MRST HO set is larger, A = 0.353 GeV when n f = 3, corresponding to a.!~ loop = 0.364 
and a? loops = 0.263 at the same scale. The shapes of the LO and NLO PDFs are also 
somewhat different which can affect results in different regions. 

A third way of defining the theoretical K factor is also possible. One could compare 
a fully NLO calculation, where the LO and NLO matrix elements are calculated with a s 
at two loops and a NLO PDF, with a LO calculation where the LO matrix elements are 
calculated with a s at one loop and a LO PDF. This procedure is possibly most favored e.g. 
in event generators where most processes are calculated with LO matrix elements only. The 
magnitude of the theoretical K factor depends on the chosen definition. Work is in progress 
to compare all three ways of evaluating the K factor for heavy quark production [ [27| |. 

Note that while the total cross section predicts the yield of heavy quark production 
over all phase space, it cannot provide much useful information on nuclear effects such a 
/^-broadening and shadowing. Any broadening will not affect the total cross section but 
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will have a strong influence on the pj distributions. Shadowing may reduce or enhance the 
nuclear cross section relative to that of the proton but the effect may be more apparent in 
some regions of phase space than others. To obtain more information on nuclear effects, 
it is thus necessary to turn to distributions. In addition, a real detector does not cover all 
phase space. The differential distributions can be tuned to a detector acceptance. For full 
details on the shadowing and broadening effects on the single quark and QQ pair distribu- 
tions, see Ref. [ ^8[]. It turns out that the pr distributions are more strongly influenced by 
any broadening effect than nuclear shadowing. Likewise, nuclear shadowing can best be 
studied combining information on rapidity and pair invariant mass distributions. The effect 
of shadowing on the dilepton continuum has been extensively studied in Ref. [ ^9|] . 

3. Relevance for Quarkonium 

To better understand quarkonium in nuclear collisions, it is necessary to have a good es- 
timate of the expected pp yields. However, there are still a number of uncertainties in 
quarkonium production in pp interactions. Two approaches have been used to describe 
quarkonium production phenomenologically — the color evaporation model (CEM) [ [$0| 
and nonrelativistic QCD (NRQCD) 

In the CEM, the QQ pair neutralizes its color by interaction with the collision-induced 
color field — "color evaporation". The Q and the Q either combine with light quarks to pro- 
duce heavy-flavored hadrons or bind with each other in a quarkonium state. The additional 
energy needed to produce heavy-flavored hadrons is obtained nonperturbatively from the 
color field in the interaction region. The yield of all quarkonium states may be only a small 
fraction of the total QQ cross section below the heavy hadron threshold, 2m//. The QQ 
cross sections we have obtained through our 'by-eye' fits of the mass and scale parameters 
have implications for quarkonium if the CEM is used to calculate the production. This is 
because different quark masses will result in more or less of the cross section below the 
2m// threshold. Since we have concentrated on bb production in these proceedings, we 
show our results for T production where there is the biggest variation in mass and scale: 
nib = 4.5 GeV is 0.85 mj while nib — 5 GeV is 0.95 me, much closer to the BB thresh- 
old. Since the NRQCD cross section is obtained independently of the parameters used to 
calculate heavy quark production, we do not discuss it here. 

At leading order, the production cross section of quarkonium state C is 



where ij — qq or gg and dij (s) is the ij — > QQ subprocess cross section. Hadronization 
is assumed not to affect the kinematics of the parent QQ pair so that only a single univer- 
sal factor, Fc, is necessary for the cross section of each state. The factor Fc depends on 
the heavy quark mass, mg, the scale /j in the strong coupling constant oc s and the parton 
densities. The factor Fc must be constant for the model to have any predictive power. The 
CEM was taken to NLO using exclusive QQ hadroproduction [ Q to obtain the energy, xf, 
and ^/--dependence of quarkonium production [ |32|, [33]]. In the CEM, gg — > g(g* — > QQ), 




(3) 
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incorporated at NLO, is similar to models of g* — > T fragmentation [ [34]]. By including 
this splitting, the CEM provides a good description of the quarkonium pj distributions at 
the Tevatron. The CEM cross sections were calculated with the MRST HO [ |l5|] parton 
densities using the same values as the NLO evaluations of the heavy quark cross sections 
shown here. The values of Fc for charmonium and bottomonium have been calculated from 
a fit to the J/\\t and T data combined with relative cross sections and branching ratios, see 
3611 for details. The combined T, T' and T" cross sections to muon pairs are 
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compared to our CEM calculations in Fig. ^. We find that the direct Y(15) production cross 
section is ss 3.4 nb at y/s = 200 GeV and varies between 126-259 nb at 5.5 TeV. 
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Fig. 4. Inclusive T production data, combined from all three S states, and compared to 
NLO CEM calculations. The solid curve employs the MRST HO distributions with nib = 
/j = 4.75 GeV, the dashed, m& = /j/2 = 4.5 GeV, the dot-dashed, m.h = 2fi = 5 GeV, and the 
dotted, GRV 98 HO with m b =/j = 4.75 GeV. 
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